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Direct evaluation of charged particle polarizabilities on the lattice is quite difficult. However, a short cut 
for charged pion polarizability - the Das, Mathur, Okubo Sum Rule - can readily be calculated using lattice 
techniques. A phenomenological model has been developed to fit the time behavior of the propagators in this 
expression. Numerical systematics are discussed and some preliminary results are presented. 



1. Introduction 

External field methods have been used previ- 
ously to measure neutral particle polarizabilities 
in the context of lattice QCD[Q. However, there 
have been no investigations of charged hadron po- 
larizabilities using the lattice techniques. This is 
unfortunate since there is now a growing body 
of experimental information on this subject. Of 
course, external field methods would be very dif- 
ficult or impossible to use in the case of charged 
hadrons because charged particles accelerate in 
an electric field. The evaluation of charged par- 
ticle polarizabilities can be done on the lattice 
by measuring a Compton scattering coefficient . 
However, the direct evaluation is actually quite 
difficult because of the many diagrams involved. 

There is a way around this problem if one is 
willing to work in the exact chiral limit. Let us 
examine a sum rule for pion polarizability and 
see how readily it can be evaluated on the lattice. 
Conventional wisdom is that charged pion polar- 
izability involves cancellation of large numbers, 
so may be difficult to simulate. 

2. DMO Sum Rule and Model 

The Das, Mathur, Okubo Smu rule for the 
charged pion polarizability, derived from current 
algebra in the chiral limit is|^: 

< > 



This employs only mesonic two-point func- 
tions, which are easily calculable on the lat- 
tice. The spectral integrals can be measured us- 
ing a combination of numerical techniques, actu- 
ally developed other physics contexts (the "vol- 
ume method" and "Fourier reinforcement" |^ ) , 
which results in a very good signal. The value of 
the term proportional to the pion squared charge 
radius, < >, and the pion decay constant, F^^, 
are taken from experiment. One can change the 
spectral integral above into a lattice space-time 
integral through the use of (the momentum trans- 
fer is purely spatial, — (0, <?)) 

r ^[p^is^)-py[s')), (2) 

where in Euclidean infinite space (Vo(a;) = 
d{x)-inu[x)) 

Koiq') = 

J d*xe-'^-- < 0\TiVoix)V^m\0 > . (3) 

(The subtracted axial propagator, A^g((f^), is 
similar to Eq.(||) but has the pion pole removed.) 
The left-hand side of Eq. (^ is formed by taking 
a numerical momentum derivative of the lattice 
data. 

There will be new sources of systematic uncer- 
tainty in the lattice evaluation of sum rules. In 
this case the above shows that one must perform 
a numerical momentum derivative. In addition. 
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the time integral can be performed in either a 
discrete or continuous sense, with attendant un- 
certainty. One can simply apply Simpson's rule 
in the discrete case. The continuous case de- 
mands some way of interpolating between the lat- 
tice time points. For this purpose, let us consider 
the standard lattice point-to-point charge density 
correlator: 

^e-^P-'<0\TiVo{x)V^m\0>. 

X 

One may show that (all masses and distances di- 
mensionless) this reduces to 

J2 e-'P-^Tr[S{x, 0)7475^^(2;, 0)7475], 



where S{x, y) is the quark propagator and the 
trace is over color and Dirac spaces. Use the free 
quark propagator (coordinate gauge; diagonal to 
this order in color space): 



1 7-2; 
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Make the replacement J^x ~* / ^"^^ ^^'^ define 



(4) 
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The next term is of order and is negligible in 
our phenomenological fits. Do the angular inte- 
grals {p=\p\): 
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r}-(6) 



d-sin(H{-^,^^,^, ^^,^^,^3^ 

Actually, what one wants in this case is the 
derivative of the above with respect to (see 
Eq.(^) above). I prefer to form the same numer- 
ical derivative in the phenomenological model as 
is necessary on the lattice. Thus, in order to com- 
pare with the lattice data, consider 



GUt.p)-GUm 
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(r2+f2)4 (^2+^2)3 
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To this order the same functional form holds for 
the axial propagator as well. 

There are two modifications to this function 
that are necessary before one compares to the 
lattice data. First, it is clear that the above ex- 
pressions have an ultraviolet infinity associated 
with the r = lower limit. This infinity can 
be controlled, as the lattice itself controls it, by 
putting in a short distance cutoff. So, replace the 
lower limit above by tq > 0, which becomes a pa- 
rameter in the fits. Call this modified function 
QQQ{t,p,ro). Second, put in a continuum thresh- 
old. So, to control the onset of excited states in 
the spectral density. One can easily show that 
the resulting function is now given by (i > un- 
derstood) 

^oo(^,-P,^o,so) = 



c 



'eo^o(t-f m,p,ro). (8) 



A third parameter, has been added in Eq.(g) 
in order to account for lattice anisotropy at small 



lattice times as in Ref . | 
the two integrals in Eq. 



L One can analytically do 
q) ; the remaining oscillat- 



ing radial integral is done numerically. One must 
also apply vector and axial vector renormalization 
constants to the lattice data before comparing to 
Eq.(^). In addition, add to Eq.(||) a function 

to fit the lowest pole mass in the vector (rho me- 
son) or axial vector (Al meson) cases. These 
mass values are fixed from previous independent 
lattice measurements, so all told there are 4 pa- 
rameters in the model: C, so and tq. 

Using Wilson fermions at /3 = 6.0 on a 16'^ x 24 
lattice, one has Pmin = f , and this value is used 
in both the numerical lattice derivative as well as 
the phenomenological model. I have not yet done 
the chiral extrapolation in quark mass, but the 
results from Eq. (|l|) when plotted as a function of 
quark mass are extremely flat. 

The fit vector correlator (unnormalized) is 
shown on a logio scale in Fig.l at k = 0.152. 
It is the area under this curve which is of interest 
here. The time integrals are very sharply peaked 
near t = 0. The numerical quality of the vector 
time data is seen to be quite good. 
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3. Summary and Remarks 



Figure 1. Vector correlator at k = 0.152. The 
dashed line is the fit from the free quark model. 



My preliminary result (at k — .154) is (using 
the lattice scale from Ref.Q) 

as = 1.7(2.1)(.8) X IQ-* fm^ 

where the first number in parenthesis is the sys- 
tematic uncertainty in the time integral, the sec- 
ond is statistical. The systematic uncertainty is 
defined to be half of the difference between the 
discrete (Simpson's rule - smaller value) and con- 
tinuous (free quark model - larger value) evalua- 
tions for as- The above is to be compared with 
the experimental result (from MARK-II data[^), 



exp 



2.2(1.6) X 10"^/m^ 



and the result from second order chiral perturba- 
tion theory 1^, 

^Mrai ^ 2.4(.5) X 10-V"^^- 

So what one gets is surprisingly good, at least 
compared to experiment. Warning: other sources 
of systematic error (momentum derivative, lattice 
scale) have not yet been evaluated. 



The sum rule method of extracting charged 
pion polarizability from lattice data has been ex- 
amined. We have seen that the systematic un- 
certainty in the time integral is larger than the 
statistical uncertainty, but still reasonable com- 
pared to experiment. 

To reduce this systematic uncertainty, the sit- 
uation here demands a fine mesh of lattice points 
in the time direction. Coarse lattices are not 
useful in the case of such sum rules. However, 
time asymmetric lattices with a small time spac- 
ing (but not too small to get into the asymptotic 
regime) could possibly be very useful. 

Pion polarizability scales like a^, and so is ex- 
tremely sensitive to the lattice scale. When the 
experimental situation is improved, it could easily 
become one of the most sensitive tests of QCD. 
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